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ABSTRACT 







A statistical analysis of 27 years of monthly averaged sea surface 
temperature (SST) and ice concentration data was conducted for 17 
locations along the annual mean position of the marginal ice zone 
spanning the North Atlantic. Anomalies (differences from monthly 
means) of both variables were observed to have spatial scales of 100s 
to 1000s of kms, temporal scales of 6 months to several years, and a 
strong regional dependence. Sea surface temperature autocorrelation 
values were in general higher than ice concentration autocorrelation 
values. Cross -correlations between the two variables were found to be 
highly significant in some regions and poor in others. The various cor- 
relation features appeared plausible with respect to understood physic 
cal processes in each region. For example, the data for the northern 
Barents and Iceland Seas showed strong cross -correlations at lags 
extending to over nine months. The steady-state cold water tempera- 
tures and relatively weak currents in these regions enhanced persis- 
tence of both SST and ice concentration, allowing them to interact. 
By contrast, the Davis Strait area, a region of strong confluent currents 
of different temperatures and limited ice persistence, showed weak 
cross-correlation values. Statistical analyses of large, homogeneous 
data sets as conducted in this study appear to be superior to current 
thermodynamic models in their potential for long-range forecasts of 
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ice concentration. 
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I. INTRODUCTION 



The ability to forecast sea ice cover in Arctic waters is of prime 
concern to many agencies. The oil industry must know the extent of 
ice to plan drilling operations and design appropriate rigs. Fishing 
fleets may be prevented from utilizing choice fishing grounds due to 
ice cover. Surface shipping must be routed to avoid the dangers of ice 
collision but retain the economy of shorter distances over great circle 
northern routes. In some regions of the Arctic, resupply by sea can 
only be conducted during minimal ice seasons. Significant volumes of 
merchant vessel traffic that operate in the northern European and 
USSR waters can only operate in relatively ice free seas. 

Sea ice introduces many complications to military operations. 
Naval vessels used in northern waters must be properly designed if 
they are-to function effectively in the unique and extreme conditions 
encountered in the region. 

The highly variable temperatures, surfaces, and material composi- 
tion of water and ice in the polar regions have a dramatic effect on 
both the air and ocean boundary layers. This in turn affects electro- 
magnetic, electro-optical and acoustic transmission and reception. 
Sea ice is a strong source of ambient acoustic noise and radar scatter. 
Detection of underwater targets by passive or active means is difficult. 
Detection of low level or surface targets by radar also becomes more of 
a problem in the ice. 
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Obviously, a reasonably precise knowledge of where the ice is and 
where it is expected to be in the future can be very useful. One 
method to approach this problem, eind indeed the focus of this work, 
is to closely examine the so-called “interannual variability” of sea ice. 
In the meteorological/oceanographic context, interannual variability 
will refer to departures from monthly and seasonal mesins defined 
from data for periods of 10 to 30 years. Longer base periods are gen- 
erally not used because climatic trends may induce changes in the 
means or “normals.” Interannual variations are the focal points of 
long-range forecasts. Long-range forecasts can be considered as pre- 
dictions of monthly or seasonally averaged departures from the clima- 
tological mean. 

In the case of atmospheric forecasts, deterministic prediction of 
daily weather fluctuations is presently limited to 5 to. 10 days. Theo- 
retical arguments indicate that the ultimate range of such predictions 
is several weeks. Evaluations of prediction skill at longer ranges is 
based on averages of forecast and observed fields over periods of 30 to 
90 days. Forecasts of 30- or 90-day average fields may well show sig- 
nificant and useful skill even though day-to-day variability is not fore- 
cast with skill. The same time scales apply to sea-ice forecasts. 

Previous studies (Lemke et al., 1980; Lebedev and Uralov, 1982; 
Walsh and Johnson, 1979) have indicated that the persistence of ice 
anomalies (departures from normal) is considerably greater than the 
persistence of atmospheric anomalies (Figure 1.1). Stated differently. 
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daily fluctuations contribute considerably more to the variance of 
atmospheric anomalies than to the variance of sea ice anomalies. 

Several other implications of Figure 1.1 deserve mention. First, 
the relatively large autocorrelations of sea ice anomalies imply that 
persistence will be a competitive control forecast. While the autocor- 
relations in Figure 1.1 are not region-specific, they imply that persis- 
tence alone can predict approximately 0.46 and 0.20 of the variance 
(autocorrelation squared) of departures from normal ice cover at lead 
times of one and three months, respectively. Persistence statistics 
will therefore be evaluated in the present work as a benchmark for 
comparison with skill statistics derived from other forecasts. Second, 
the sea ice autocorrelations in Figure 1.1 are comparable to those of 
sea surface temperature (SST), which has long been used as an input 
to long-range weather forecasts for mid-latitudes. The similarity of 
time scales suggests that SSTs are associated with ice concentration 
anomalies. If the two variables are coupled, then the skill of predic- 
tion of one variable may be enhanced by considering the distribution of 
the other. 

Variable sea surface temperatures and ocean currents have been 
cited as likely sources of the “unexplained” ice variance in studies of 
atmospheric forcing of sea ice (Haupt and Kant, 1976). Schell (1970) 
and Rogers and van Lx)on (1979) have compared large-scale SST and 
sea ice variations in selected regions and selected years. Considerable 
effort has also been expended in evaluating the coupling between 
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large-scale SST and atmospheric anomalies over the North Pacific 
(e.g., Namias, 1976; Elsberry and Raney, 1978; Broccoli and Hamack, 
1981) and the North Atlantic (e.g., Ratcliffe and Murray, 1970; 
Rowntree, 1976; Haworth, 1978). Because these studies have indi- 
cated the existence of statistically significant associations between 
SST anomalies and the atmospheric circulation, the link between SST 
anomalies and the high-latitude sea ice distribution must be viewed in 
the context of Figure 1.2. The “direct association,” in which high 
SSTs enhance melting and/or retard freezing, is to be distinguished 
from the “indirect association,” in which the sea ice distribution is 
influenced by SST through an SST-induced effect on the atmospheric 
circulation. The interpretation of the SST- ice coupling in observa- 
tional data analysis is complicated by the possibility that both mecha- 
nisms may act concurrently. Indeed, while the present study will 
focus on the predictive applications of the “direct” association 
between SST and sea ice in the North Atlantic, the “indirect” associa- 
tion is the basis of Soviet claims concerning seasonal predictability of 
North Atlantic sea ice (Chapter IV). 

The present work is a data-based investigation of SST-derived 
predictability of sea ice in the North Atlantic. The study is motivated 
by three major considerations: 

1. SST represents potentially valuable but as yet largely unexploited 
input to long-range sea ice prediction 

2. Previous studies addressing the high-latitude portions of the 
Pacific and Atlantic Oceans have achieved considerably less 
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diagnostic success in the Atlantic than in the Pacific (e.g., Walsh 
and Sater, 1981); and 

3. The recent release of the Comprehensive Ocean Atmosphere 
Data Set (GOADS) provides the most solid foundation to date for 
an analysis of SST/ice concentration associations. However, the 
potential for high latitude forecasting is still unclear as GOADS 
has yet to be utdized for this purpose. 

The following points will serve as objectives for this study: 

1. Autocorrelations of sea ice and SST anomalies will be systemati- 
cally evaluated in order to place the SST-derived ice pre- 
dictability into a framework of practical utility 

2. Lag relationships between SST and sea ice concentration 
anomalies will be evaluated in order to determine which variable 
(SST or ice concentration) may serve as a useful predictor for 
the other in high-latitude regions. 

3. Soviet studies (e.g., Lebedev and Uralov, 1982) contain claims of 
high levels of skill in seasonal forecasts of summer ice coverage 
in the Greenland Sea, Since an SST input contributes to this 
skill, the Soviet scheme will be tested and used as one bench- 
mark for the evaluation of the sea ice concentration-SST associ- 
ations obtained here. 

4. In view of (3), particular attention will be given to the spatial and 
seasonal generality of the SST/ice concentration relationships 
obtained for regions such as the Greenland Sea. Specifically, the 
temporal and spatial scales of North Atlantic sea ice concentra- 
tion and SST fluctuations will be evaluated. This should allow 
assessment of regional and seasonal dependencies of SST- 
derived forecast sldll. 

The statistical correlations describing the SST/ice coupling will 
be diagnosed by isolating the contribution of individual years to the 
correlations, and by then examining the seasonal evolution of the 
anomalies within these years. 

The methodology chosen for use in this study is primarily statisti- 
cal rather than model-oriented. The choice of this strategy was made 
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because numerical models have not yet proven effective for long-range 
ice forecasting, even when run uncoupled to variable SSTs. The use of 
numerical models for operational forecasting of ice conditions by the 
United States Navy is presently limited to the 0-6 day range, and this 
model (see Preller, 1985) uses prescribed ocean currents and heat 
fluxes. State-of-the-art coupled ice-ocean models of the North Atlantic 
are constrained by either the use of prescribed North Atlantic fluxes 
into the Arctic (Semtner, 1987) or by a damping of the ocean variables 
to a prescribed climatology (Hibler and Bryan, 1987). Even when 
these constraints are eventually eliminated, the ice/ocean models will 
require prognostic forcing fields (at least from the atmosphere) for 
applications to monthly or seasonal ice forecasting. The relatively low 
levels of skill shown by atmospheric forecasts at monthly or seasonal 
ranges implies that numerical ice model forecasts for these ranges will 
be severely limited by atmospheric input. It is therefore likely that 
the present reliance on statistical techniques in monthly and seasonal 
ice forecasting (e.g.. Stringer et al., 1984) will continue through the 
foreseeable future. This was the main consideration in the decision to 
use statistical techniques for this work. 

Chapter II contains a description of the ocean characteristics of 
North Atlantic high latitudes, with an emphasis on the variables rele- 
vant to seasonal forecasts of the ice cover. A description of the data 
sets used and manipulations applied to them follows in Chapter III. 
Chapter IV contains an evaluation of a Soviet study of Greenland Sea 
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Chapter IV contains an evaluation of a Soviet study of Greenland Sea 
ice predictability. The results provide one of the benchmarks for 
comparison to the results obtained in the present work. Chapter V 
presents the results of this thesis and is followed by a discussion in 
Chapter VI. Finally, the conclusion and recommendations for future 
work are included in Chapter VII. 
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AUTOCORRELATION 




Figure 1.1 Compilation of autocorrelations of SST, sea ice, snow and 
700 mb pressure amomalies from various authors. Auto- 
correlations were obtained over large areas of the north 
Pacific (SST and 700 mb) North American continent 
(snow) and north Atlantic (ice). Data were averaged over 
seversd months, [from Davis (1976), Lemke et al, (1980), 
Walsh (1982) and Namias and Bom (1970)1. 
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Figure 1.2 Schematic representation of the coupling mechanisms between sea surface temperature 
(SST) and sea ice extent [from Walsh and Sater (1981)]. 



n. ARCTIC OCEAN CHARACTERISTICS 



A feature which contributes to the uniqueness of the oceanogra- 
phy of the Arctic region is sea ice. Coachman and Aagaard (1974) 
state: 

The general oceanographic consequences of a perennial or seasonal 

ice cover are: 

1. The water temperature of the near-surface layer in the presence 
of ice is always maintained close to the freezing point for its 
salinity by the change of phase process. 

2. Salt is excluded from the ice to a varying extent, but the water 
under the ice is always enriched in salt by any ice growth. The 
dependence of water density on temperature and salinity is such 
that close to the freezing point density is almost solely a function 
of salinity. Therefore, ice formation can increase the density 
locally and some vertical convection may result. 

3. In the transfer of momentum from the atmosphere to the ocean, 
the wind must act on the sea through the intermediary of the ice. 

The focus of this work is on the relationship between SST and ice 
concentration. Points 1 and 2 listed above indicate two mechanisms 
which serve to make the relationship an important one. The third 
point, especially in the resultant form of ocean currents, is also 
relevant. 

A closer look at these three factors (ice, water and currents) is 
considered useful to put this study in perspective. Much of the infor- 
mation contained in the following sections A through C was extracted 
from Coachman and Aagaard (1974); Sater et al. (1971); and CIA 
(1978). 
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A. ICE 



The Arctic Ocean ice pack is confined by a nearly continuous 
boundary of land. The associated constraint on equatorward transport 
is a major reason why ice which forms here survives longer and devel- 
ops into more complex forms than ice found in southern polar regions. 
An annual net heat loss and stratification of the underlying water also 
contribute to Arctic ice longevity. 

Ice formation in open water starts in the autumn. As days grow 
shorter and nights longer, the amount of solar insolation decreases. 
Since radiation from the earth remains approximately constant, the 
energy budget changes from a net gain in summer through equilibrium 
to a net loss in the fall. The relatively warm mixed ocean surface layer 
of summer is cooled until reaching its freezing point and ice crystal 
formation commences. 

Ice initially forms around the boundaries of the polar ice pack and 
over the shallow protected waters of high-latitude coastlines. The 
ice-covered areas continue to expand until they merge, forming an 
ice-locked Arctic Ocean from October to June. The ice cover starts 
from a summer minimum of approximately 5.2 million square km, 
then more than doubles in areal extent to a maximum of 11.7 million 
square km by the end of the ice season (CIA 1978, p. 12). 

The ice cover grows continuously until spring. Once the days 
start to lengthen, the amount of solar radiation increases until the 
surface energy budget is agciin positive. Snow cover will reflect up to 
90 percent of the sun’s incident radiation, which slows the initial 
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heating stages. However, once the air temperature reaches the snow’s 
melting point, the albedo rapidly drops to approximately 40 percent. 
This results in a period of rapid snow melting. 

Continued heating initiates melting of the underlying ice, causing 
cracks and flaws to develop. The surface melt water drains through 
the cracks, further eroding them until the ice breaks up into floes. 
Eventually, the ocean surface layer heats up, the floes gradually dis- 
solve, and the cycle is complete. 

The extremes of ice cover in the period 1953-1977 are shown in 
Figure 2.1. Minima and maxima for both summer and winter are pre- 
sented. Obviously, sea ice cover has a large degree of variability both 
inter-seasonally and inter-annually. Ice cover is something of a mis- 
nomer as the sea surface is rarely covered by an unbroken expanse of 
ice. Even in the very thick ice regions of the winter polar pack, 
infrared measurements indicate that up to 10 percent of the area of 
the ocean is either open water or thin ice from recently refrozen leads 
(Sater et al., 1971, p. 41). 

Multi-year ice (ice which has survived a summer’s melting) com- 
prises the majority of the polar pack, which averages 3 m in depth 
(Bourke and Garret, 1987). First-year ice rarely grows to a thickness 
greater than 2 m. However, the depth of the ice in any location is 
largely dependent on the external forces of wind and currents. These 
cause the ice to converge and diverge. When a region of ice con- 
verges, it buckles, folds, and overlaps, forming a rugged terrain and 
areas of considerably thick ice. For example, the Beaufort Gyre’s 
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anticyclonic flow causes ice to converge along the north coast of 
Ellesmere Island and Greenland. The number of ridges in this region 
is well above the average for the Arctic pack (Weeks, 1978). The 
mean ice thickness here is of the order of 6 to 8 m (Bourke and 
Garrett, 1987). In contrast, the ice pack east of Spitzbergen 
(Svalbard) is not confined by land and is free to diverge. Herer the 
average ice thickness is significantly less, averaging approximately 
2 m. 

B. WATER MASSES 

The waters of the Arctic Seas are often described on the basis of 
temperature and salinity. As such, they are comprised of three main 
water masses: Arctic Surface Water, Intermediate or Atlantic Water, 
and Deep or Bottom Water (Coachman and Aagaard, 1974). Arctic 
surface water is generally limited in depth to 200 m. It has the most 
variable characteristics and can be modified by the weather, the sea- 
son, and/or the physical environment. Temperatures in this layer vary 
from -1®C to over 2®C. The salinity may be uniform to approximately 
50 m, below which a sharp halocline increases the salinity to about 
34.5 %o at the bottom of the layer. The variety of conditions in the 
surface water is evident in Figure 2.2, where vertical profiles of tem- 
perature and salinity from a variety of Arctic basins are plotted. 
Coachman and Aagaard (1974, p. 9) state that 

The most important processes conditioning and modifying the sur- 
face layer are: 



25 



1. Addition of mass (fresh water) from the land, primarily from the 
large Siberian rivers; 

2. Additions of fresh water locally through melting of ice; 

3. Heat gain through absorption of solar radiation in non-ice- 
covered areas during summer; 

4. Concentration of salt and hence increase of density of surface 
water, through freezing of ice; 

5. Heat loss to the atmosphere through any open water surface, 
including leads in the central Arctic pack ice; and 

6. Inflow and subsequent mixing of Atlantic and Pacific waters. 
Processes 1, 2, and 3 occur only from June to September and lead 

to a decrease in water density. These buoyant waters form a surface 
cap which absorbs radiation and warms. Therefore, in summer, ice- 
free regions tend to have warmer and less saline surface layers. In 
areas where the ice does not recede, surface temperatures remain 
near freezing as incoming energy is used to melt the ice but surface 
layer salinities are reduced due to ice melting. 

Processes 4 and 5 have the greatest impact in winter. In some 
areas, such as the shelf waters, conditions are such that the water 
becomes dense enough to penetrate into the intermediate Atlantic 
layer (Aagaard et al., 1985). However, in general, the strong pycno- 
cline at the base of the surface layer prevents mixing due to surface 
density changes from penetrating below 200 m (Coachman and 
Aagaard, 1974). 

As a consequence of the halocline-derived mixing barrier, pro- 
cess 6 does not have a great influence on the surface water variability. 
However, regions of high surface salinities (33 to 34.5 %o) are found in 
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the Greenland and Labrador Seas and East Baffin Bay. These high 
salinities are due to the advection of North Atlantic surface water into 
the Arctic via warm surface currents from the south. One further con- 
sequence of the halocline barrier is that, in winter, when the surface 
layer cools below freezing, the relatively warm Atlantic layer is insu- 
lated from the ice cover. This severely limits any vertical heat flux. 
However, cases have been observed along the continental slopes where 
the Atlantic layer is forced to shallower-than-normal depths. Vigor- 
ous surface mixing can then break through the halocline, and vertical 
heat flux can occur (Coachman and Aagaard, 1974). This results in the 
ice melting from the bottom or not forming at all, forming open poly- 
nas in the ice cover. 

The second water mass is called the Atlantic Layer. It is the 
result of the influx of warm and salty (35.0 to 35.1 %o) North Atlantic 
water flowing into the Arctic basin through Fram Strait and the 
Barents~Sea (Weigel, 1987). It extends from 200 m to 900 m with 
temperatures above 0® C. A temperature maximum of approximately 
0.45° C, observed throughout the Arctic Basin, occurs between 300 m 
and 500 m, dependent on location. The salinity gradually decreases to 
approximately 34.9 %o in the Arctic Ocean and Greenland Sea and 
approximately 34.6 %o in Baffin Bay. 

The final water mass, the bottom water, has temperatures below 
0° C. The salinity is nearly constant from the bottom of the Atlantic 
layer to the ocean floor. The intermediate Atlantic water and deep 
water are advected into and out of the Arctic seas from adjacent areas. 
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principally through Fram Strait in the North Atlantic. The difference 
between intermediate and deep water is simply defined as where the 
temperature is above 0® C and below 0° C, respectively. Both water 
masses are nearly isohaline and therefore of nearly uniform potential 
density. 

A temperature difference does exist between deep waters of^ the 
Canadian and Eurasian basins. The deep Canadian basin averages 
approximately -0.45® C while the Eurasian basin temperature is 
approximately -0.79® C (Aagaard et al., 1985). The two temperatures 
are kept separate by the Lomonosov Ridge, which acts as a sill 
between the two disparate water masses. 

Observations of the vertical structure at various stations have been 
taken over long time periods and in many different locations. A 
remarkable similarity in the profiles has led to the conclusion that the 
Arctic basins are in a long-term d 5 niamic steady-state condition 
(Coachman and Barnes, 1961). It is further noted that observed dis- 
tributions of Arctic water properties are a result of continuing pro- 
cesses within the basins. Therefore, surface water T-S profiles 
indicate the local modifying processes, while T-S profiles for depths 
below 200 m indicate the common origin of the water. 

C. CURRENTS 

The surface circulation of the Arctic basin has been derived from 
satellite observations and the plotting of ice island, buoy and floe sta- 
tion movement, as well as ship’s tracks. The circulation of the Arctic 
waters is due both to water density differences and wind forcing. 
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Although large anomalies in the flow (compared to the long-term 
mean currents) are often observed, the long-term mean currents have 
been reasonably well documented in the north Atlantic seas (Krauss, 
1986). A chart covering the majority of the area of interest in this 
study is included in Figure 2.3 which indicates the primary ocean 
currents. . 

Circulation and ice movement patterns in the peripheral Atlantic 
seas of the Arctic Ocean are dominated by two major currents systems, 
one warm, the other cold (Krauss, 1986). The warm North Atlantic 
Drift (NAD) Current splits off the Gulf Stream extension south of 
Newfoundland and heads northeast. It divides at the Mid-Atlantic 
ridge at approximately 51° N. One portion flows north to Iceland, 
where it becomes the Irminger Current. The Irminger Current itself 
then splits, with one branch taking water eastward around the north 
coast of Iceland, where it becomes the North Icelandic Current (NIC) 
and then the East Icelandic Current (EIC). The other branch flows 
westward, then to the south along the east Greenland slope, eventually 
rounding Kap Farvel and heading into the Davis Strait. The main por- 
tion of the NAD passes through the Scotland-Faroes gap and up along 
the Norwegian coast, where it is known as the Norwegian Atlantic 
Current. At the northern tip of Norway, the current splits again. One 
portion continues along the north shore of the USSR to Novaya Zemlya 
as the North Cape Current, while the other portion, called the West 
Spitsbergen Current (WSC), heads north into the Arctic Ocean just 
west of Svalbard. 
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The Transpolar Drift Stream is a major cold current originating 
under the polar ice cap and exiting the Arctic Ocean through Fram 
Strait, between Svalbard and Greenland. It passes down the east coast 
of Greenland as the East Greenland Current (EGC) and splits Just 
north of Iceland (Aagaard and Coachman, 1968). One portion contin- 
ues down the Greenland coast, rounds its southern tip, and heads 
north up the western coast of Greenland as the West Greenland 
Current. This portion merges with water flowing through the various 
passages of the Canadian Archipelago to form the cold Labrador 
Current. The Labrador Current follows the east coast of Canada from 
Baffin Island to Newfoundland. The second of portion of the split 
north of Iceland, known as the Jan Mayen Current, flows 
southeastward until it merges with the EIC and then finally with the 
northeastward-flowing Norwegian Atlantic Current. The southward 
EGC, eastward Jan Mayen Current, and north and westward flowing 
WSC combine to form a cyclonic gyre in the Greenland Sea. Filaments 
of cold water from the Transpolar Drift Stream also flow south from 
the Arctic Ocean in a broad, slow pattern from Svalbard to Svemaya 
Zemlya. This flow essentially limits the northward warming influences 
of the Norwegian Atlantic Current. 

Figure 2. 1 clearly shows the impact of the major current patterns 
on the mean ice coverage. The cold currents support ice growth and 
enhance the spatial extent of sea ice, while the warm currents melt 
the ice or preclude its formation. The southerly extent of the ice edge 
between regions of dissimilar temperatures can vary as much as 30 
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degrees of latitude. For example, the cold Labrador Current can often 
support ice cover as far south as Newfoundland. However, the North 
Cape of Norway, which is under the influence of the warm Norwegian 
Atlantic Current, rarely experiences sea ice. An extreme example is a 
region northwest of Svalbard, where open water can normally be found 
in wintCT as far north as 80® N. This is over 800 km farther north 
than any other open water in this season (Sater et al., 1971). 

The currents also explain many other features of ice coverage and 
movement. Sea ice off the east coast of Greenland is found to originate 
primarily in the Arctic basin and as a result can be quite thick. The 
East Greenland Current continually carries a wide belt of Arctic pack 
ice southward through Fram Strait. As the ice season develops, the ice 
edge migrates farther and farther south down the east Greenland 
coast. Extreme years will have the ice edge as far south as Kap Farvel 
or as far north as 70® N. Although sea ice exists north of 70® at all 
times oTihe year, the belt narrows and has more open water in the 
late summer (Sater et al., 1971). 

Strong contrasts in ice coverage are also observed within the 
Baffin Bay-Davis Strait region. The warm WGC follows the bathymetric 
contours of the southern and western Greenland continental shelf, 
keeping the southwest coast of Greenland ice free during most win- 
ters. However, the cold Labrador Current, as noted before, contributes 
to the large concentrations of ice found along Canada’s east coast well 
south of Baffin Island. This current also carries the majority of 
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icebergs southward through the Grand Banks fishing zones and the 
Hibernia oil fields, and into the north Atlantic sea lanes (CIA, 1978). 

D. SEA/ICE FORCING 

The primary driving forces for sea ice growth and movement are 
the air stress resulting from wind-induced surface drag and the water 
stress resulting from ocean currents. Variations in wind forcing are 
determined by the passage of s 3 moptic scale weather systems having 
spatial scales of 100 to 1000 km and a duration of several days. 
Responses of the sea ice to mesoscale features of the atmospheric and 
oceanic forcing produce eddies and other fluctuations with length 
scales of 10 km and time scales of hours to days. Inertial oscillations 
of the drift of sea ice have similar periods. Longer term fluctuations, 
with time scales of weeks to seasons, represent responses to atmo- 
spheric and ocean forcing integrated over equivalent time scales. 

When examining fluctuations var 5 dng on time scales of a season to 
several decades, the dominant signal in sea ice variability is clearly the 
annual cycle (Parkinson et al., 1987). Figure 2.1 contains 25-year 
envelopes of the positions of the Arctic ice edge at the end of February 
and August, the approximate times of maximum and minimum ice 
extent. At most longitudes, the seasonal change from summer to 
winter ice edge positions is considerably greater than the 25-year 
range of extremes for a particular month. The annual cycle is primar- 
ily a response to the seasonally varying insolation and may therefore be 
attributed to thermod 5 mamic rather than dynamic forcing. 
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Figure 2.1 Ice extremes for the period 1953-1977 (from Walsh and Johnson (1979)). 
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2 Vertical profiles of temperature and salinity from various 
northern high-latitude basins, [from Coachman and 
Aagaard (1974)]. 
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Figure 2.3 Schematic picture of the North Atlantic circulation as 
derived from drift experiments. Vertical shading indicates 
flow associated with the North Atlantic Current and hori- 
zontal shading indicates association with the subtropical 
g 3 rre. Area between A and A' is the shifting boundary 
between the two. Warm currents are marked by solid 
arrows and cold currents by dashed arrows, [from Krauss 
(1986)]. 
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in. DATA 



The statistical analysis in this study were based on two data sets. 
A modified COADS (Comprehensive Ocean-Atmosphere Data Set) data 
base provided monthly averaged sea surface temperatures, and the 
SEIC (sea ice concentration) data base provided sea ice concentrations 
(in tenths) determined for the end of each month. 

A. COADS 

COADS contains approximately 70 million reports from ships of 
opportunity, ocean weather ships, buoys, and bathythermographs. The 
data have been assimilated, sorted, edited, and summarized statisti- 
cally for each month during the years 1854 to 1979 in 2® latitude X 2® 
longitude boxes over the entire world’s oceans (Slutz et al., 1985). 
Data sources are noted in Table I. 

COADS is considered to be the most complete data set of envi- 
ronmental parameters now available for the ocean/atmosphere bound- 
ary. However, as described below, Slutz et al. (1985) suggested that 
careful attention must be paid to the sometimes serious, and usually 
poorly understood, limitations inherent in any data set of this type. 

Potential error sources for the COADS data set include: 

1. Instruments, observation and coding methods, navigation accu- 
racy, ship construction, and data density have aU undergone his- 
torical changes. Many of these changes were unrecorded in the 
data sets from which COADS was derived. 

2. SST measurement. Temperatures measured by engine intake 
have been shown by Ramage (1984) to be approximately 0.5® C 
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warmer than temperatures measured by bucket, but Slutz et al. 
(1985) indicate that, in a data set in which both bucket and 
intake temperatures are used indiscriminately, there is currently 
no way to effectively differentiate between the two and make 
appropriate corrections. Fortunately, inhomogeneities 
attributable to this problem were not severe during the time 
period examined by this study, 

3. Diurnal effects. The data sets identify how much daylight had 
occurred prior to each observations. However, effects discussed 
by Ramage (1984), such as insolation, heating from the ship,- and 
cloud cover biases, are not taken into account. 

4. Significant errors can occur at every stage of observation, 
recording,transmission, and processing. Outlier filters were 
applied, but these can only remove the most obvious mistakes. 

5. Duplication. Data files often contain numerous duplicate entries 
of one observation. This obviously biases the data. Over 25 per- 
cent of the original data were rejected because of duplication. 
However, other duplications undoubtedly still remain. 

The region of interest for this thesis was contained in latitudes 
28° North to 88° North. Preliminary review indicated that the com- 
pression of the longitude lines as they converge at the pole reduces 
the area of the 2° X 2° boxes to the point where the number of data 
reports per box become too sparse for reasonable analysis. As a result, 
the GOADS data were averaged into 4° longitude X 4° latitude boxes. 
The average sea surface temperature in each box was then assigned to 
a grid point at the center of the 4° box. This resulted in a monthly 
grid of 90 X 15 = 1350 points. 



B. SEIC 

The SEIC data set containing monthly ice concentration was 
compiled by Walsh and Johnson (1979) from existing sea ice distribu- 
tion observations. The original data sources are noted in Table II. The 
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data set was recently updated with the digitized ice data for 1972- 
1984 from the United States Navy/NOAA Joint Ice Center (Gross, 
1986). 

These data sets were compiled by applying observations to the 
1648 data points contained in a 58 x 80-point grid. The spacing was 
1° latitude (approximately 110 km). The grid covered those portions 
of the Arctic seas where ice was observed during any month. The 110 
km spacing was chosen to permit resolution of year-to-year fluctua- 
tion. Grid point values contain the fraction (tenths) of area covered by 
ice within each box. 

The data set contains values for each month from 1953 to 1984. 
1953 was chosen as the first year since data records for several 
regions of the Arctic seas became essentially continuous from that 
period on. 

Walsh and Johnson (1979) noted that SEIC, like the GOADS, had 
limitations. Briefly, these included: 

1. Tlie data set covers a time period both before and after satellite 
coverage. Different observational methods cause nonuniformify in 
the data interpretation. 

2. Imprecise concentration classifications. Ice conditions can vary 
considerably over relatively small areas. Ice-observing agencies 
such as the British Meteorological Office and the United States 
Naval Oceanographic Office tend to group the concentrations into 
categories such as “very open” (1/10 to 3/10) and “open” (4/10 
to 6/10). When the source charts were digitized, the center 
value of these classification groups was used. 

3. Overlapping data. When more than one source covered the same 
area, discrepancies were sometimes observed. In these cases, 
the mean was the value digitized and applied to the grid point. 
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4. Missing data. Some points in the early years had no data, partic- 
ularly the Siberian Sector during the 1950s. In such cases, grid 
values were either linearly interpolated between prior and fol- 
lowing months or a 25-year (1953-1977) monthly mean was 
used. These estimated values were flagged in the digital storage 
mode. 

C. DATA SET MANIPULATION 

Correlations between SST and ice concentration required _data 
covering the same years. The time period of January 1953 to 
December 1979, a total of 27 years or 324 months, was selected for 
both variables as each was available over this time frame. 

The 58 X 80 ice concentration grid was considered the simplest 
to work with. Because the SST data field was reasonably dense, a sim- 
ple bi-linear objective analysis was used to shift the SST values to the 
58 X 80 grid. The four adjacent SST values were weighted as a linear 
function of distance. Data points on land were compensated by apply- 
ing a zero weighting factor to land values but maintaining the total 
“correctILweightings at 1.0. 

Statistical analysis of every grid point was both impractical and 
unnecessary. The prime concern was making predictions in areas of 
the marginal ice zone (MIZ). Therefore, a procedure was designed to 
select designated points of interest (POI) in that region. 

Correlations of SST with regions of near ice absence (Q^io) or total 
ice coverage generally do not produce good results as the ice 

concentration does not vary much at these extreme conditions. The 
most useful correlations are obtained where large changes in ice con- 
centration occur. The entire 324 months of ice concentration values 
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were averaged in order to form a mean grid for the entire sample 
period. Since the annual cycle contained the largest ice concentration 
variance, the 5/10 average concentration contour on this mean grid 
was considered the most probable location for these large variance 
values. A total of 17 well-spaced POIs were chosen along the 5/10 
contour from Hudson’s Bay to the northern tip of Norway. The points 
were selected to cover a wide variety of current, latitude, water mass 
type, and data density regimes. Subsequent analysis revised these 
POIs to permit better usage of available data and more consistent 
spacing north of Norway. The final POIs, plotted on the ice 
concentration grid, are shown in Figure 3.1. Ice concentration values 
at each POI were then averaged with values from 24 surrounding 
points. The 24 points were those which had positions two or less grid 
steps in both directions from the POI. This 25-point average provided 
some smoothing to the data. 

Monthly averages for each POI over the 27 years were calculated 
for both ice concentration and SST. These monthly means were then 
subtracted from each mean monthly POI value to produce a monthly 
anomaly field. This process removed the annual cycle leaving monthly 
anomalies of SST (STDMA) and ice concentration (ICDMA). The 
removal of the annual cycle severely reduces the ice/SST correlations 
which, in “raw” form, are dominated by the annual cycle. The use of 
departures from the monthly mean focused concentration on the 
interannual fluctuations of the two variables. 
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Initial correlation runs were conducted using no time averaging, 
three-month centered averaging, and five-month centered averaging. 
Increasing the time averaging reduced the noise, allowing better 
determination of the correlation signal. Consequently, a five-month 
centered averaging was chosen as the final smoothing routine. A 
precedent for the use of five-month smoothing was established by 
Climate Analysis Center (NOAA/NWS/NMC). Its monthly climate diag- 
nostics bulletin (NOAA, 1987) contains a number of products (e.g., the 
Southern Oscillation Index) where five-month averaging has been used 
to good effect for reducing contamination of signals by short-term 
noise. A summary of the data set manipulations is included in Figure 
3.2. 
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